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Abstract 

The renormalization of the most general dimension-six four-fermion operators without power 
subtractions is studied at one loop in lattice perturbation theory using overlap fermions. As 
expected, operators with different chirality do not mix among themselves and parity-conserving 
and parity-violating multiplets renormalize in the same way. The renormalization constants of 
unimproved and improved operators are also the same. These mixing factors are necessary to 
determine physical matrix elements relevant to many phenomenological applications of weak 
interactions. The most important are the K'^ — K^ and B^ — B° mixings in the Standard Model 
and beyond, the A/ = 1/2 rule and e'/e. 



1 Introduction 

Since the original proposals of using lattice QCD to study hadronic weak decays Q , substantial 
theoretical and numerical progress has been made. In the most popular lattice regularizations, i.e. 
Wilson and staggered fermions, the main theoretical aspects of the renormalization of composite 
four-fermion operators are fully understood ||2|, ^]. The calculations of the matrix elements relevant 
for the K^-~K^ and B^-B^ mixings have reached a level of accuracy which is unmatched by any 
other approach [Q, ||]; increasing precision has also been achieved in determining the electro- weak 
penguin amplitudes necessary for the prediction of the CP-violation parameter e'/e p|-[12]. Finally, 



matrix elements of AS* = 2 operators, relevant to the study of FCNC effects in several extensions of 
the Standard Model (supersymmetry, left-right symmetric models, multi-Higgs models etc.), have 
also been computed |l^]-|jl3|. Some of these lattice predictions have been fundamental in constrain- 



ing the parameters of the CKM matrix in the Standard Model (SM) |14| , 15| and beyond |16]. 

Nevertheless, some fundamental phenomena, such as the A/ = 1/2 rule in non-leptonic kaon 
decays, and the value of e'/e, which measures the direct CP violation in kaon decays, are far from 
being understood. In the Standard Model these effects can be explained only if the non-perturbative 
physics gives contributions to matrix elements definitively larger than their factorized values ||l5|| . 
Therefore a non-perturbative determination of the relevant matrix elements is crucial for predict- 
ing these quantities. Lattice QCD is the only method which can address these problems from first 
principles. Techniques have been developed for both Wilson and staggered fermions, but these 
methods have not yet produced useful results |^, 17|. 



To compute non-leptonic weak matrix elements it is essential to construct renormalized oper- 
ators in definite chiral representations. In the Wilson or SW-Clover lattice regularizations, bare 
operators do not have a definite chiral behavior due to the presence of the symmetry breaking term 
in the action. Renormalized operators with the correct chiral properties are recovered as linear 
combinations of operators with different chirality. The Standard Model O ^^ operator is the 
most popular example: its matrix element between pseudoscalar states should vanish as M^ in the 
limit Mtt -^ 0, while the matrix elements of the wrong-chirality operators which mix with it are 
expected to go to a constant (and in fact in the kaon mass region they are 2 to 10 times larger 
than the SM one). Therefore the correct chiral behavior of O "^ is recovered only if the finite 
mixing constants are known with high precision. This has been a long-standing problem and has 



only been solved using Non-Perturbative (NP) renormalization techniques y, |l^, 18 1. The situation 
becomes even worse when the lack of chiral symmetry complicates mixings with lower dimensional 
operators. This is one of the obstacles to a reliable computation of the matrix elements relevant 
for A/ = 1/2 and e'/e. On the other hand, with staggered fermions chiral symmetry is preserved, 
but solving the doubling problem and defining operators with the correct flavor and spin quantum 
numbers is far from trivial. 



Only recently has it been understood [|19|]-|23] that chiral and flavor symmetries can be preserved 



simultaneously on the lattice, without fermion doubling, if the fermionic operator D satisfies the 
Ginsparg- Wilson Relation (GWR) @: 

75Z? + Z?75 = -D-i^D . (1) 

P 



The GWR implies an exact symmetry of the fermion action at non zero-lattice spacing, which may 
be regarded as a lattice form of the standard chiral rotation |23|. Nevertheless, it is important to 



stress that locality, the absence of doubler modes and the correct classical continuum limit are not 
guaranteed by the GWR in eq. (jTI). Indeed, there exist lattice fermion actions which satisfy the 



GWR but which do not meet the above requirements |25|. A breakthrough in this field was achieved 
through the domain- wall formulation of lattice fermions |26] and by Neuberger through the overlap 
formulation [^I|. He found a solution of the GWR which satisfies all the above requirements and 
is locafi 13 : 

Dn = ^(l+X^l= 
a V VXfX 

X = Dw--P, (2) 

a 

where D\y is the Wilson-Dirac operator and < p < 2r (see below). A further remarkable property 
of Ginsparg- Wilson fermions is the absence of 0{a) discretization errors in the action and therefore 
in the spectrum of the theory. Nevertheless, the local fermionic operators have to be improved to 
remove the 0{a) effects in matrix elements. This step is greatly simplified by the Ginsparg- Wilson 



relation, as it allows the construction of 0(a)-improved operators to all orders in g^ |28|, which 
renormalize with the same renormalization factors of the unimproved ones. On the other hand, the 
complicated form of the Neuberger operator in eq. ^ renders its numerical implementation quite 
demanding for the present generation of computers. However, some progress has been achieved |2^ 
and Monte Carlo simulations seem to be already feasible, at least in the quenched approximation. 

Once the action satisfies the GWR and all the properties described above, the quark masses 
renormalize only multiplicatively and mixings of operators with different chirality are forbidden 
even at finite cut-off [^]. The most general set of dimension-six four-fermion operators without 
power subtractions is 

Or^rs = V^ir^V'2 • V53rU4 (3) 

Or^rs = V^ir^V'4 • v53r'BV2 . 

The main result of this paper is the evaluation of the mixing pattern of these operators at one 
loop in perturbation theory in the overlap lattice regularization: we compute their renormalization 
constants and we show how the Neuberger action greatly simplifies the mixing among them. 

The renormalization constants are necessary for extracting physical matrix elements from nu- 
merical simulations. The matrix elements of the operators studied in this work are relevant in 
many phenomenological applications of weak interactions. They are necessary for predicting the 
K^ — K^ and B^ — B^ mixing amplitudes in the Standard Model and in several extensions of it (su- 
persymmetry, left-right symmetric models, multi-Higgs models etc.). By using chiral perturbation 
theory, they can also be related to the A/ = 3/2 matrix elements relevant for the prediction of e'/e. 
They are also necessary for estimating the B^ — B^ width difference and the 0{l/'m^) corrections 

^From eq. (H) it is clear that the overlap operator is not ultra-local. The Neuberger kernel satisfies a more general 
definition of locality, i.e. it is exponentially suppressed at large distances with a decay rate proportional to 1/a. 



in inclusive b-hadron decay rates |^^. The renormalization constants we have computed are also 
relevant for the study of the AI = 1/2 rule and e'/e on the lattice. The operators involved in these 
computations can mix with same- as well as lower-dimension operators. The key observation in 



these cases |31] is that the lower-dimension operators do not change the anomalous dimensions of 
the original operators. Therefore the mixings among dimension-six operators can be obtained from 
those of the operators studied in this paper. 

Perturbative computations with the overlap-Dirac operator are much more cumbersome than 
for Wilson fermions, due to the complicated structure of the vertices in the action (see Appendix 
^. For this reason it is very useful to use an intermediate renormalization scheme to separate as 
much as possible the computations performed using lattice perturbation theory (where, for example, 
the Fierz rearrangements work) from the continuum perturbation theory which is much simpler. 
The RI renormalization scheme proposed in ||3^, ^ is the optimum choice: we first renormalize 
lattice operators in the RI scheme using only lattice perturbation theory and then the matching 
with one of the MS schemes is done in continuum perturbation theory. Lattice and continuum 
regularizations are used independently taking full advantage of their properties. 

In Ref. [p^ j, renormalization constants of local bilinear fermion operators have been computed 
in the MS renormalization scheme using the lattice overlap regularization. We computed these 
renormalization constants in the RI scheme, matched them in MS and verified that our results are 



in agreement with those of Ref. |34]. Then we computed the renormalization constants of the four- 
fermion operators we are interested in following two independent ways: the direct computation, like 
in Ref. |35], and one which uses the Fierz and charge conjugation rearrangements to compute the 



four- fermion Green's functions in terms of the bilinear ones ||36| , g . The agreement of the results is 
also a non-trivial check that the numerical integrals with the Neuberger propagators and vertices 
are well estimated. 

A further determination of the renormalization constants considered in this paper can be 



obtained using numerical NP methods such as those in Refs. |33, B9l, applied to Neuberger fermions. 
However, perturbative estimates are very useful because often they are very good approximations 
and they furnish a consistency check of the NP methods. Moreover, for Neuberger fermions the 
perturbative computations can remain the only determinations of the renormalization constants 
for some time. 

The paper is organized as follows: in section y we define the overlap fermion action used; in 
section |3| we discuss the bilinear renormalization constants; in section ^ we address the mixing of 
the four-fermion operators and report the results for the renormalization constants; in section |5| we 
state our conclusions. 



2 Basic definitions 

The QCD lattice regularization we consider for massless fermions is described by the action 
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Sl = -E 1 - 7:Tr [C/p + 4] +J2J2Mx)DN{x,y)My), (4) 



i=i x,y 



where, in standard notation, Up is the Wilson plaquette and qq the bare couphng constant. D^ is 
the Neuberger-Dirac operator defined in eq. (|2|), where the massless Wilson operator Dv/ is defined 
as 

^w^ = ^7/.(V^ + Vp-^aV;V^ (5) 

and V^ and V* are the forward and backward lattice derivatives, i.e. 

V^V'(a;) = - Uij,{x)il){x + aji) -ii{x) 



a 



V!l^(x) = - iljix) - UUx - afi)'ilj{x - afi) . (6) 



The range of the Wilson parameter is < r < 1 and U^{x) is the lattice gauge link. Eq. (|^) implies 
a continuous symmetry of the fermion action in (^, which may interpreted as a lattice form of 
chiral symmetry [p^]: 

5^ = 75 M - —aDTvj ■0, (5^ = V5 f 1 - — oDatJ 75. (7) 

The corresponding flavor non-singlet chiral transformations are defined including a color group 
generator in eq. (^. The generalization to massive fermions is simple [^0|: in eq. (Q) Dat has to be 
replaced by 

Dn — >\l- ^arnM Dn + ml (8) 

where rrig is the bare physical quark mass of flavor i. The Feynman rules of the action deflned in 
Eq. (Q) are given in Appendix]^ and are in agreement with those used by |^, |4l|, ^]. Throughout 
this paper we will use only mass independent renormalization schemes (RI and MS) and therefore 
all our computations are performed with massless quarks. 

3 Renormalization of Bilinear Operators 

In this section we set our notation for bilinear operators and we report the results we have obtained 
for their renormalization constants in the RI and MS schemes. 
A generic non-singlet quark bilinear is deflned as 

Or{x) = Mx)TM^) (9) 

where the flavors ipi — 1^2 are different and E is a generic Dirac matrix. The GWR ensures that 
no lattice artifacts of 0{a) are present in the action and therefore also not in the spectrum of 
the theory. However, matrix elements of operators are still affected by 0{a) discretization errors 
that have to be removed by improving the fermionic operators. In Refs. |2^ it is shown that, for 
massless quarks, the improved bilinear operator is given by 

1 1 

Of = V^(l - ^a^Jv) r (1 - —aDN)^ . (10) 



It can be proven [34, 28] that the renormahzation constants of Of are the same as those of the 
corresponding Or- Therefore ah the results obtained in this section for local bilinear operators also 
remain valid for the corresponding improved operators. 
The renormalized bilinear operators are defined by 

Or(/i) = Zr{afi)Or{a) (11) 

and the chiral symmetry in eq. (^) imposes the constraints 

Za = Zv 

Zp = Zs (12) 

on the renormalization constant^. 

The quark propagator in momentum space is denoted by S{p) (for the conventions adopted 
throughout the paper see Appendix y). The two-point fermionic Green's function of a bilinear 
inserted at the origin {x = 0) is 

Gr(xi,X2) = {M^i)Or{0)Mx2)) , (13) 

its Fourier transform with equal external momenta is 

Grip) = J dxidx2e-*P(^^-^2)Gr(xi, X2) , (14) 

and the corresponding amputated correlation function is defined as 

Arip) = S-\p)Grip)S-\p) . (15) 

The renormalized quark propagator is 

S{p,fi) = Z^{afi)S{p,a) , (16) 

1/2 

where ZJ is the quark field renormalization constant, and the renormalized Green's functions are 

Gr(p,/i) = Z^{fj,a)Zr{fia)Gr{p,a) , (17) 

A^{p,^^) = Z^H/xa)Zr(/xa)Ar(p,a). (18) 

The anomalous dimensions of composite operators are defined as 

7r(a.) = -2Z,-V^^r = E7?(^) , (19) 



^If the conserved axial and vector currents corresponding to the chiral symmetry in eq. (pn were used, then 
Zy = Za = 1 (l|. 



where Ug = fl'^/47r is the strong couphng constant. If a renormahzation scheme which preserves the 
chiral Ward Identities (WI) is chosen (such as MS or RI), then 



IP 



7V = , 

75 = -7m , 



(20) 



(0)c 



where 7m, is the anomalous dimension of the quark mass. Moreover the 7p s are renormahzation- 
scheme independent and gauge invariant and their values are reported in Table |l|. 

The anomalous dimensions are properties of the renormalized operators: they are independent 
of the regularization scheme adopted and depend only on the renormahzation conditions imposed 
(the renormahzation scheme) and the renormalized coupling constant used. For a given operator, 
defining the renormahzation conditions is equivalent to fixing 7r(as), i.e. 7r(Q;s) uniquely deter- 
mines the scheme in a regularization invariant way. Therefore, 7r can be computed in a simpler 
regularization, like Naive Dimensional Regularization (NDR), provided that the renormahzation 
conditions are the same. The evolution of the renormalized operators is determined by solving the 
Renormahzation Group Equations and, in a given scheme, is 



Or(^) = #^Or(/xo) 



cr(/^o) 



(21) 



where cr(/x) at Next-to-Leading Order (NLO) 



(0), 



cr 



(/.) = (a.(^))-^V(2A,) 



1 + 



Q^(m) (tv A7r 



47r 



2/3o 



2/3o^ 



(22) 



is scheme dependent. In eq. (^) /3o and /3i are the first two coefficients of the /? function of QCD. 
For 7p in RI or MS see for example [|^]. The Renormahzation Group Invariant (RGI) operators 



can be defined as |43]-|45| 



0?«i 



cr(/^) 



Or(/i) 



Zf«i(a)Or(a) 



where 



Zr\a) 



Zrjpa) 
cr(^) 



(23) 



(24) 



Op*^^ is independent of the scheme, scale and gauge chosen to renormalize the operators, while 
Zp (a) depends only on the regularization and not on the scheme, scale and gauge. 

As discussed in the introduction, it is very helpful to separate computations which can be 
performed using only lattice perturbation theory from those done in the continuum. In this regard 
the RI renormahzation scheme proposed in ||3^ , |33| is the optimum choice: it preserves all the 
relevant symmetries (chirality and switch, see below) and it can also be easily implemented non- 
perturbatively [^]. The matching with a "continuum" renormahzation scheme, for example MS, 
remains necessary because almost all the Wilson coefficients are computed in the continuum, but it 
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S 


V 


T 


Continuum PT 


^(0) 





-6 





2 


^MS,RI 





-4 










P 


Lattice PT 




0.2 


-235.80762 


1.31942 


1.52122 


1.58848 




0.3 


-150.61868 


1.89625 


1.52277 


1.39828 




0.4 


-108.19798 


2.38060 


1.52448 


1.23911 




0.5 


-82.86081 


2.80522 


1.52637 


1.10009 




0.6 


-66.05227 


3.18782 


1.52845 


0.97532 




0.7 


-54.10921 


3.53927 


1.53074 


0.86124 




0.8 


-45.20179 


3.86686 


1.53329 


0.75543 




0.9 


-38.31447 


4.17577 


1.53611 


0.65622 


B 


1.0 


-32.83862 


4.46989 


1.53924 


0.56236 




1.1 


-28.38734 


4.75224 


1.54274 


0.47290 




1.2 


-24.70304 


5.02527 


1.54665 


0.38711 




1.3 


-21.60760 


5.29104 


1.55105 


0.30438 




1.4 


-18.97397 


5.55135 


1.55601 


0.22422 




1.5 


-16.70910 


5.80783 


1.56163 


0.14623 




1.6 


-14.74330 


6.06201 


1.56804 


0.07005 




1.7 


-13.02336 


6.31544 


1.57541 


-0.00460 




1.8 


-11.50798 


6.56970 


1.58394 


-0.07798 



Table 1: Perturbative values in the continuum and on the lattice in the Landau gauge, where 
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(0) 
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^^' /Cp while the Bs contain the proper vertex contribution only. 



can be done using continuum perturbation theory only. Therefore lattice and continuum regular- 
izations are used independently taking full advantage of their properties. In the following we will 
indicate the renormalization scheme adopted with a superscript in the renormalization quantities, 
i.e. Zy and Zp^^ for RI and MS schemes respectively. 

Using the Feynman rules defined in Appendix ^, we computed the self energy of the quark 
propagator and the amputated Green's functions of bilinear operators between off-shell quark states 
at one loop in perturbation theory in a generic covariant gauge. The full expressions we obtained 
are reported in Appendix |D| and they can be used to impose any renormalization condition in any 
covariant gauge. The RI scheme is imposed by taking the trace of the amputated Green's functions 
in the Landau gauge (a = 0) ||3^. The renormalization constant of the quark field is then defined 
as 






(25) 



while the renormaUzation constants of the bihnear operators are given by 



(Zf(/ia))-iZ«-i(/xa)Tr PrAr(p,a 



(26) 



where the trace is over both color and spin indices and Pp are suitable projectors on the tree-level 
operators defined in Appendix |y. In general, these renormalization conditions differ from the stan- 
dard MOM perturbative prescriptions (see for example |^6[), where the renormalization constants 
are extracted only from terms proportional to the tree-level matrix elements of the operators under 
consideration. The differences are evident from the expressions of the quark propagator and the 
amputated Green's functions reported in Appendix ^. In the standard procedure, finite terms 
such as ^Pfi/p'^ in the vector current are considered as part of the matrix elements of the opera- 
tors, whereas with the projectors they give additional finite contributions to the renormalization 
constants^. These terms are the same on the lattice and in the continuum and they cancel when 
one computes the difference between the continuum and the lattice renormalization constants, as 
in the standard perturbative procedure. 

Using eqs. (25) and (p6|) and the results reported in Appendix O we obtained 



7-ni 



Z^' 



167r2 



Cf 



^-li^^- 



-(0) 



-(0) 

7r 



log(/ua)2 + B^ 



log(^a)2 + Bt + B^ 



(27) 



where the 7p and Bs are reported in Table ||. 

Once the operators have been renormalized in RI, the matching with another given scheme x 
is obtained with 

Z^{„a) = 4{i.)Zf^\a) = ^Zf\„a) , (28) 

Cp (^/ij 



where the Cp^(//) are reported in |43] for the bilinears, defined with the same conventions adopted 
in this paper. 



To compare our results with |p4|, we report explicitly the matching coefficients for the MS 
scheme. The chiral WI imply Z^^ = Z^, Z^^ = Z^ and 



7RI 7MS 



A straightforward computation in NDR gives 



1 . 



(29) 



Zf\a^) 






Z^\afi), 



(30) 



^Note that in the Landau gauge these terms are absent both in the bihnears and in the four-fermion operators. 



where the A^^'^^s are reported in Table |l|. In the MS scheme, the renormahzation constants 
of gauge-invariant operators are independent of the gauge. Therefore, the dependence on the 
gauge and external states of the RI scheme cancels the corresponding dependence of the matching 
coefficients. Our results in MS are in agreement with those of Ref. [p4| . 

4 Four-fermion operators 

In this section we introduce the most general set of dimension-six four-fermion operators we are 
interested in, we analyze their mixing pattern exploiting only the symmetries of the underlying 
theory and finally we compute the necessary renormahzation constants at 1-loop in perturbation 
theory. A very exhaustive non-perturbative analysis for Wilson fermions has been done in Ref. p!q| . 
We will proceed on the same lines with the advantage of having the additional chiral symmetry in 
eq. (0) which forbids mixing among operators which belong to multiplets with different chirality 

The generic set of four-fermion operators we are interested in is 

Or^rs = ^ir!4V'2 • V53rU4 (31) 

Or^Ts = V^ir!4V'4 • v53rU2 , 

where the fiavors tpi — ipi are all different and the T^ are generic Dirac matrices with / representing 
the contracted Lorenz indices (if any). The 20 operators in (31) form a complete basis of four- 



fermion operators. Operators with the color indices contracted in a different way can be expressed 
as linear combinations of the ones in eq. (^) using the color Fierz identity (^) in Appendix ^. 

The 0{a) discretization errors in on-shell four-fermion matrix elements can be removed by 
using the improved operator 

^(l _ ^uDn) Ta (l - ^aDjv)^ • ^(l - y'^Dn) Tb (l - ^aDjv)^- (32) 



This can be proven along the lines used in Refs. |28] for the case of two-quark operators. Here 
we just sketch the argument, which is valid for any reasonable action which satisfies the GWR. 
Starting from the Neuberger operator D]\f, one can define the associated operator |25] 



Kn = {i- ^«^w) ^Dn, (33) 

which has the same chiral properties as the continuum Dirac operator, i.e. {-f^iv,75} = 0. The 
propagator K^ is free of 0{a) corrections and, although otherwise not well-behaved and so not 
useful in practice, it turns out to be very useful for the construction of operators that are improved. 
In fact, the four-fermion correlation function 

1 /. 1 „ \„ /. 1 „ \ 1 1 /. 1 „ \„ /. 1 „ \ 1 \ 



_^l _ -aD,)T,(l - -aD,)—^—(l - ^^aD,)T,(l - ^^aD,,—, 



is automatically improved. In the last line we have used eq. ( |33[ ) to re-express the Green's function 
in terms of the Neuberger propagator, and from this we can read off the improved operator. Thus, 
the Ginsparg- Wilson relation highly simplifies the improvement of the four-fermion operators at all 
orders in perturbation theory. Moreover, the renormalization factors for corresponding improved 
and unimproved operators are the same. This happens because in 1-loop amplitudes a Djsi factor 
can combine with a quark propagator but, since it has an a in front and (contrary to the Wilson 
case) there is no 1/a factor in the propagator, as additive mass renormalization is forbidden by 
chiral symmetry, its contribution to the renormalization factors is zero |34|. 

The symmetries relevant for studying the mixing of the operators that we consider are Parity 
("P), the Switch (5) transformation ^2 '^ ^4 and chiral symmetry in eq. ^. They allow splitting 
the original basis into smaller independent multiplets. For parity-conserving ("P = 1) operators it 
is useful to define the following bases 



and for the parity-violating i(P 



where 



Qf 


= 


O^v + Oaa 


Qt 


= 


Otv - OtA 


Qt 


= 


0^ n^ 


Qt 


= 


^SS + '-^PP 


Qf 


= 




= -1) 






Qf 


=^ 


OtA + O^y 


Qt 


= 


OtA - Oty 


Qt 


= 


-Ogp + Opg 


Qt 


= 


Ogp + Opg 


Qt 


-- 


0±~, 




TT ' 


Ot.r, = 


1 
2 


Or,r, ± 0^^ 



(35) 



(36) 



(37) 



and T = cr^u'^^. Since the lattice preserves parity, the two sets of operators (^) and (|3q ) renormalize 
independently. Using continuous chiral transformations it is easy to show that, at variance with 



Wilson fermions, they renormalize with the same renormalization matrices [18|. We have explicitly 
verified this property at one loop in perturbation theory. In the following we will then consider 
only the parity-conserving operators in eqs. (^). Among them, the five Qf are left invariant under 
a switch transformation [S = 1), while the Q~ change sign {S = —1). Therefore the two sets 
renormalize independently as 



Qf{^i) = Z^i^ia)Qfia) 



(38) 



where Z^ are the renormalization constant matrices. 

The operators of the bases (|35|) do not belong, in general, to irreducible representations of 



10 



the chiral group. Nevertheless chiral symmetry imposes many constraints on the matrix of the 
renormahzation constants when a scheme which preserves it together with the switch symmetry is 
chosen. The most general mixing matrix under these constraints reads 



(39) 



Note that for i,j = 2,3, Z^ and Z-- are related. Therefore for Neuberger fermions there are only 
14 independent renormahzation constants, which would instead become 64 if the Wilson action 



/7± 











\ 





Z22 


±^23 












I 


±Zs2 





^33 







^54 




Z55J 



were usedQ. The particular structure of the matrix in eq. (3£) also simplifies the implementation 
of non-perturbative techniques |3^, ^ . 

Analogously to the bilinears, denoting by xi, 2:3 and X2,X4 the coordinates of the outgoing and 
incoming quarks respectively, the four-point Green's functions are defined as 



Gf{xi,X2,X3,X4) = {'lpl{xi)lp2ix2)Qf{0)lp3ix3)lp4.{x4,)), 



(40) 



where (• • •) denotes the vacuum expectation value. Note that Gf depends implicitly on the four 
color and Dirac indices carried by the external fermion fields. The Fourier transform of the Green's 
function (|^) at equal external momenta p is defined as 



Gf{p) = / dxidx2dX3dx4e-^P(^i+^3-^2-^4)Gf (xi,X2,X3,X4) . 



The corresponding amputated correlation functions is defined as 

Af{p) = S~Hp)S-\p)Gt{p)S~Hp)S~\p) . 

From the above definitions we find for the renormalized Green's functions 

Gfip) = Z^Zpfip), 
Afip) = Z;^'Z^Af{p). 

The anomalous dimension matrices are defined as 



(41) 



(42) 



7^(as) 



d 



-2(Z±)-V-^Z±=E7^« 
"*" i=0 






j+i 



(43) 



(44) 



*In the Wilson case the renormahzation matrices of parity-conserving and parity- violating operators are different 
[iq ]. Moreover all the entries of the analogous of the matrices in eq. (p9|) are independent and non-zero. This 
results in 25 X 2 + 14 = 64 independent constants. It is interesting to note that using suitable Ward Identities, some 
parity-conserving matrix elements can be related to the parity- violating ones |46|. 
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and, if the renormalization scheme preserves chiral symmetry, they have the same structure as in 
eq. (|39|). At first order in Og we obtain 

(45) 



7ii — 37s 


7ii 


4_(o) 


7S^ = -^# 


723 - 


-27^ 


4^ = 


(0) 
733 = 


hP 


+(0) 23# + 9# 
744 - 12 


744 


_ 4l7f - 97^ 
12 


-(0) , -(0) 

.+(o)_ -7s +7t 

745 - 12 


-(0) 
745 


_5(#-#) 
12 


..(0) mP-iP) 
754 - 4 


754 


-7™ +7?' 

4 


+(0) -9#+41# 
55 1 


755 


_ 94°) + 237?) 

1 o 



which agree with |3^, ^, where the anomalous dimension matrices at two loops in the same bases 
can also be foundQ. 

The renormalization group evolution of the four-fermion operators, in a scheme which preserves 
chiral and the switch symmetry, is 

Q±(/i) = w^{is)[id^if,o)]-'Q^if^o) , (46) 

where in this case w (fi) are matrices which depend only on the anomalous dimensions and as- 



Their expression at NLO in a generic scheme can be found in |37, 38, 12 1. Analogously to the 



bilinears, the RGI operators can be defined as |12| 



[Q±]RGi = [^±(^)]-iQ±(^) ^ [Z±(a)]^GiQ±(a) , (47) 

where 

[Z^{a)fGl ^ [y;±(^)]-i^±(^a) . (48) 

As for the scheme, scale and gauge dependences, these definitions have the same properties as for 
bilinears. 

We have computed the amputated Green's functions of the four-fermion operators at one loop 
in two different ways. In the first one we use the Fierz and charge conjugation rearrangements to 
connect the proper vertices of the four-fermion operators to the bilinear ones. In the second, the 
calculation has been performed computing the four-quark diagrams directly without any rearrange- 
ment of quark legs. To carry out the analytic calculations we have used a set of routines written 
in the symbolic manipulation language FORM, which are a generalization to the overlap case of 



the ones used in Refs. [35|. Many of the features of those calculations are also present here. We 



^The operators Qg used in this paper and defined in equation (BSl) correspond to 1/4 of those defined in M, Kq]. 
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Figure 1: An example of Fierz and charge conjugation rearrangements for Feynman graphs of 
four-fermion operators. 




CHARGE 
CONJUGATION 




FIERZ 
REARRANGEMENT 




have numerically integrated the analytic FORM outputs for some values of the p parameter and 
we have successfully checked the results obtained with the first method. This also turns out to be 
a non-trivial check of the numerical integrations relevant in the Neuberger regularization. The full 
expressions for the one-loop amputated Green's functions that we have obtained in a generic covari- 
ant gauge are reported in Appendix 0, and they can be used to impose any set of renormalization 
conditions. 



Once the results have been rotated in the basis (p3), the RI renormalization conditions are 



imposed by taking the trace the amputated Green's functions in the Landau gauge [|18| 

6ik {i,k = !,■■■ ,5) , 



T^Pf^k 



(49) 



where the projectors Pj are defined in Appendix |y. Using the wave function renormalization Z^ 
defined in eq. (|2^) and the definition in (E3) we obtain 



(50) 





[Z^f' = 1 - 


9l 

16^2 


"7(0) 

2 log(/^«) 


+ B^ + 2CfB^ 


where 


B^ = 


(Bf, 


: 


I 




B22 
^Bs2 







±-823 

B33 






\ 



-B44 ^45 

-^M B^^ J 


5 


and the B-- are 

'J 












^" 3 


'^By) 


^u 


_ A{Bs + By) 

3 




R ^^^ 
B22- 3 + 


3Bv 


B23 


= -2{Bs - By) 




B^ 


2 = 


B33 


8Bs 
3 
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(51) 



(52) 



+ _ 23Bs + 9Bt _ 4.1Bs - 9Bt 

_ -Bs + ^r „_ _ HBs - Bt) 
^45- 12 ^45" 12 

^+ 5(Bg - Bt) „_ -^5 + Bt 

^54 - ^ ^54 - ^ 

, _ -Ws + AIBt „_ _ 9^s + 2Wt 

where the i?r's are reported in Table |l|. The relations above are valid in the Landau gauge for any 
action which preserves chiral symmetry, where Ba = By and Bp = Bs- 

The Wilson coefficients and the four-fermion operators of the weak effective Hamiltonian are 
often defined in one of the MS schemes. Unfortunately the definition of the MS scheme for the 
four-fermion operators is not unique even if we consider the Naive Dimensional Regularization 
only. This is a source of confusion in the literature, where quite often one encounters comparisons 
of matrix elements that have been computed in different schemes. Incidentally, we note that the 



MS-NDR scheme used in the lattice calculation of Ref. [Q] differs from the one used in |37], which 
moreover is not the one adopted in |^^. In some cases, differences between various schemes may 
be numerically large. On the other hand the matrix elements of the four-fermion operators in a 
given scheme are useful only if they are matched with the Wilson coefficients computed in the same 
scheme. From the perturbative computation of the Wilson coefficients in a given scheme x, it is 
straightforward to extract the evolution operator [{v^{iJ,)]'^ at the NLO. Then the corresponding 
renormalization constants are given by 

[Z^if^a]]^ = [^±(M)]^[Z±(a)]^«i = [w^{^,)nw^''\f,)]~'[Z^{f,a)f' , (53) 

where [w (fi)] can be found i n [p7| , [l2| , p8[ . This is equivalent to using only RGI Wilson coefficients 
and operators as proposed in |l2]. 



5 Conclusions 

The matrix elements of dimension-six four-fermion operators without power subtractions are the 
primary ingredients for studying many interesting phenomenological problems in weak interactions, 
among which the most important are the predictions of the K^ — K^ and B^ — B^ mixings in the 
Standard Model and beyond. The renormalization constants of these operators are also relevant 
for studying the A/ = 1/2 rule and computing the CP-violation parameter e'/e. 

In this paper we have studied in detail the mixing of these operators using the overlap lat- 
tice regularization proposed by Neuberger and computed their renormalization factors at one loop 
in perturbation theory. The computations were done in two independent ways: in one case, us- 
ing Fierz and charge conjugation rearrangements, the four-fermion Green's functions are given in 
terms of the one loop results for the bilinear operators; in the other case the Feynman diagrams 
are computed directly using FORM codes. 

We have shown that operators belonging to different chiral representations do not mix among 
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themselves. We have exphcitly verified that, as expected, the mixing matrices for the parity- 
violating and parity-conserving operators are the same. Furthermore, the improvement of the 
matrix elements is highly simplified by the Ginsparg- Wilson relation and the renormalization con- 
stants are the same as the unimproved case. In the Wilson formulation, whether improved or not, 
the construction of renormalized operators requires subtractions of operators with wrong chirality. 
These become very severe in the case of power divergences, i.e. A/ = 1/2, e'/e, etc. We believe 
that the overlap regularization represents a very promising approach for solving these long-standing 
problems. 
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A Feynman Rules 

In this appendix we report the Feynman rules used in the computations. Some abbreviations of 
functions occurring on the lattice are useful: 

a z g 

1 ' 



k^ = - sm{k^a) k^ = J2kl. (54) 



« 



The gluonic propagator in a generic covariant gauge is defined as 

G^^Ak) = (5^, + (a - l)k^K/P) /P , (55) 

where a is the gauge parameter. 

Let us define some useful matrices in Dirac space: 

(56) 



Xoip) -- 


= h^lP^^ 


_ p 

a 


ar ,9 




Xi^{p + p') = 


= -igo 


fP^^+p'^ \ . . fP^,+p'^, \ 
7m cos y ^~~^"j ~ *^ '^^ V 2^ / 




^2m^(P+/) = 


- "4^ 


rcosy ^ "^aj Z7^sm(^ ^ '^aj 
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The fermionic propagator of the overlap action is 



^S'w^^ 



-^IpPp 



uj{p) +b{p) 



+ 1 



a Xlip)+u;{p) 
2p uj{p) + b{p) 



with 



b{p) = y^ ~ ~ ' 



uj{p) 



A^)i 



/p-2+(^p2_^ 



(57) 

(58) 
(59) 



For < p < 2r the propagator S}^ (p) exhibits a massless pole only when p^ = 0, and there are no 
doublers. 



The vertices of the overlap Dirac operator are 



and 

v!.i\p,p') 



uj[p')+u}{p) I "^ ' uj{p') ^' uj{p) J 






(60) 



+ 



2{a;(p') + uj{p)}{uj{p) + w(p + k)}{uj{p + k) + uj{p')} 
X {Xo(/)^iV (.P + P+ k) Xi, {2p + k)+ Xi^ {p + p' + k) Xl{p + k)Xi,{2p + k) 
+ Xi^ {p + p' + k) Xl {2p + k)Xo (p) 



(61) 



io{p') + (jj{p) + uj{p + k) 
uj{p')lo{p)lo{p + A;) 



Xo{p')xl {p + p' + k) Xoip + k)xl {2p + k) Xoip) 



B Notations and Fierz Transformations 



In this appendix we fix our notations for the color and spin matrices. We also report the for- 
niute used for Fierz transformations in color and Dirac space. 

The Gell-Mann group generators of the SU{Nc) Lie algebra are denoted by t , A = 1, . . . , N^ — 1. 
They are Hermitian, traceless Nc x Nc matrices and are normalized according to 



^ ^ 2 



They satisfy the commutation relations 



[t^^t^]=ifABC^C 



(62) 



(63) 
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where summation over repeated indices is implied. The structure constants / are completely 

antisymmetric and real. With these conventions the completeness relation for the t-generators 
reads 



A=l 
and using the above formulas we get 



\6(i )cd 



-Z [ Sadhc 



-TT^ab^cd 



N'~l 



A.A\ 



Er^ 



Jab 



CfK 



A=l 

E / 

C',D=1 



ACD fBCD 



'f 



Ca6 



AB 



(64) 

(65) 
(66) 



Nc. 



(67) 
(68) 
(69) 



with Cf = (iVc - l)/2iVc and Ca 

The complete basis of 16 Euclidean Dirac 4x4 matrices is denoted by 

r = {1, 7m' ^/^i^' 7m75, 75} = {5', V, T, A, P}, 
where ■~fi are the usual Euclidean Dirac matrices in four dimensions and 

The charge conjugation matrix satisfies 

and in our basis C = 7072- 

Repeated T matrices imply summation of their Lorenz indices (if any); for example VV = 
J2riln '^ 7/i5 ^^ = Jl u<u '^ fj.v ^X" (T^jy, where the sum is over the 6 independent a^i, matrices only. 
With these conventions the T matrices are normalized as 

SS =1, VV = 4, TT = -6, AA = -4, PP = 1 

where summation over Dirac indices is understood. 

The Fierz transformation of the Dirac indices of a four-fermion operator is defined as 

r r EE r„/3 ® r^5 ^ [r r]^^ = r^s ® r^p . (7i) 

The Euclidean Fierz-transformed Dirac tensor products [F F] ^ can be re-expressed as a linear 
combination of the complete set of the original tensor products F F as follows : 



(70) 



[V vf^ 
[T Tf"" 
[A Af'' 

V [P0Pf^ J 

The overall minus sign is due to the anti-commutativity of the Fermi fields. 



/ 



V 



1 


1 


-1 


-1 


4 


-2 





-2 


6 





-2 





4 


-2 





-2 


1 


-1 


-1 


1 



1 \ 

-4 
-6 

4 

1/ 



/ S0S \ 

V(E)V 
T(g)T 
A^A 

V p^p J 



(72) 
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C Projectors for Bilinears and Four-Fermion Operators 

In this appendix we define the projectors used to impose the RI renorniaUzation conditions for the 
bihnear and four-fermion operators. 

For bilinears the projectors are defined as 

Ps = — / Pp = — 75 

b 12 12 ' 

Pt = -^cj^u, (73) 

where the identity color matrix is not shown. 

If we define 

Prr = r r , (74) 



its trace on a generic four-fermion amputated Green's function is defined as |18] 



TV PrrArr(p) = T^p <S) T „, p, Krr{p)'^^pt^\ (75) 

where the color and spinor indices are explicitly reported and the trace is taken over spin and color. 
The projectors P^ for the parity conserving operators Q^ are |18| 



^' = +64iv^(kTT)(^-^^-) 



P^ ^ +7r77T?^-^^i^vy " ^ Aa) ± :jT71^777f^^ (^55 - ^PP 



1 

64(iV|"- 1) '^ " '^ ^ ^^' ^ 2,2N^{N^ - 1) 

Pf = ^ 32iV.(^.^-l) (^^^ - P-) + 16(iVf3Ty (P- - P-) (^6) 

^ ~ 32A^,(A^2 _ 1) ^^'ss + fpp) + ^2N,{N^ - 1)^^ 
P± = ^ i (P.o + Ppp) + (2A^cTl) p 

They are obtained from the tree-level amputated Green's functions A^ imposing the following 
orthogonality conditions: 

Tr Pf aJ(°) = 5,fc {i,k = l,...,5). (77) 



The analogous projectors for the parity- violating operators can be found in |18| 



18 



D One-loop results 



The one-loop expression for the quark propagator in a generic covariant gauge is 



s-^ = i:^ \l + 



167r2 



Cf 



aL + B^ + ad, 



i-tp 



(78) 



and the amputated Green's functions of the bihnear operators in a generic covariant gauge are 



As = S{l + 



al 



■Cf 



(3 + a)L + Bs-a {d.^ - 1)] } 



167r2 

Kv = v{l + -^CF\-aL + Bv-adA]-2-fj^CFa^ 



167r^ 
9l 



p. 



At = T[l + ^CF[{l-a)L + BT-a{d^ + l)]} 

2 2 

Aa = A[l + -^CF[-aL + BA-ad^]}-2-^CF 

(3 + a)L + Bp-a {d^ -!)]}> 



(79) 



A, 



167r2 



p2 



with L = log(pa)2, and the B^s are reported in Table ||. The terms proportional to the gauge 
parameter a can be given in simple form as functions of d^ = 7^; — Fq — 1 = —4.792010 (where ^e 
is the Euler constant and Fq = 4.3692252 • • •). They depend on the gluonic action, and they are 
independent of the p parameter |42], as can be seen using gauge Ward Identities. 

The four-fermion operators that we have considered in the "direct" calculation has the form 



OtaTb 

Ope pc 






(80) 



The one-loop expressions for the amputated Green's functions of the parity-conserving operators 
in a general covariant gauge are 



Asc 



Os^s^ + -^l [ - ^( - 2(3 + a)^ + ci - 2a{d^ - 1)) - N,(^aL + C3 + a(d^ + ^))] Os-5- 



iV. -1//1 



1 



PaPf3 



+ "yv^s ((9-^ + C2)0tT - -aapa «> (Tpf^ 



A^c - 4 / / 1 



N, VV2 



(^(^-L + C2J0t':t'= - -aapa_t^ ® Oppt 



p2 



Ass = Oss + -^l^^^^j^{-2i3 + a)L + ci-2aid^-l))Oss 
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Ape 



Opcpc -\- 



167r2| 



2iV. 



2(3 + a)L + ci - 2a{d^ - 1) ) - NJaL + C3 + a[d^ + -] ] Opcp 



iV. -1//1 






+ j^2 iio + ^2 - 2«JOtT + 2"^pa (^ ^P/3 



A 



pp 






Av<= 



ycyc 



01/=^= + 



167r2| 



2iV. 



2aL + C4 — 2ad^ 1 — iVc((- + alL + C6 + ad^ 



Ol/c 



v^yc 



+ { — + NAaf- :,^^- 



+- 






iV| 



-L + C5 -ajO^A + a 



p2 



iV^//3 N |^75^^^J^t-^75 



A- 



vv 



Ovv + 



9l fiVc^-1 
167r2| 2Nc 



2aL + C4 — 2a(i^ 1 Oyy — 4a 



i)^i) 



+4(^(|-L + C5 - ajOA<:A<= + a 



1^75^^ ^ ^75^^ 
p2 



Aa- 



Oa'^A'^ + 



167r2| 



1 



+ (^ + iVc]a 



2iV, 
J^75^^ ^ 1^75^^ 

p2 



2aL + C4 — 2a;d^ ] — NAi- + a\L + cq + ad. 



Oa 



cA-: 



+nv^VV2 + ^5 - aj Ovv + a ^^ 



+- 



iV. -4//3 



Nr 



-L + C5 — a ) Oy^y^ + a 



p2 



^AA = Oaa + 



al 



2aL + C4 — 2ad.j, ) O^A — 4a 






1^75 ® 1^75 



+4((-L + C5 - a)Ov<=y<= + a ■ 



pz 



(81) 
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At^c 



Ot^t^ + iI^I 



2N,, 



2(1 - a)L + CT- 2a{d^ + 1)) - iVc((2 + a)L + cg + a(d.^ - ij))\Ot-T' 



Ni-l 



Ni 



+ '^2 (^-^ + ^8 - 2") ('^^^ + ^^^) + "j^ (3^ + C8 - -a) (O5-S- + Op.pc) 

Ktt = Ott + yI^ I ^^^^(2(1 -a)L + C7-2a(dv, + l))OTr + 4(3L + 
The Cj constants, which take the same values also for the mixings of the parity-violating operators 



C8 - 2") {Os'^S" + Opcpc 





p = 0.4 


p = 1.0 


p=1.7 


Cl 


4.761193 


8.939783 


12.630887 


C2 


-0.142685 


-0.488442 


-0.790006 


C3 


-1.524484 


-1.539241 


-1.575409 


C4 


3.048969 


3.078482 


3.150818 


C5 


-0.428056 


-1.465325 


-2.370017 


C6 


-1.952540 


-3.004566 


-3.945426 


C7 


2.478227 


1.124716 


-0.009205 


C8 


-0.856112 


-2.930650 


-4.740034 


Cg 


-2.095225 


-3.493007 


-4.735432 



Table 2: Values of the finite constants in Landau gauge. 



(as we have explicitly verified), are reported in Table ^. They correspond to the values of the finite 
parts in Landau gauge and, in terms of the finite parts of the bilinears, are: 



Cl = 2Bs 
a = 2Bv 
c-r = 2Bt 



C2 
C5 

cs 



Bt - Bs 



By - Bs 



C3 
C6 



Bs + 3Bt + 4Bv 



Bs + By 



-^(Bs-Bt) C9 = -^i3Bs + BT). 



(82) 



The coefficients of the logarithms can also be re-expressed in terms of of the anomalous dimensions 
of the bilinears. 



References 



[1] N. Cabibbo, G. Martinelh, R. Petronzio, Nucl. Phys. B244 (1984) 381; 

R. C. Brower, M. B. Gavela, R. Gupta, G. Maturana, Phys. Rev. Lett. 53 (1984) 1318; 
C. Bernard, Argonne 1984, Proceedings of Gauge Theory on a Lattice, p. 85. 

[2] M. Bochicchio et al., Nucl. Phys. B262 (1985) 331. 

21 



[3] S. Sharpe et al., Nucl. Phys. B286 (1987) 253. 

[4] G. Kilcup, S. Sharpe, R. Gupta, A. Patel, Phys. Rev. Lett. 64 (1990) 25. 

[5] S. R. Sharpe, Nucl. Phys. B (Proc. Suppl.) 53 (1997) 181; 
L. Lehouch, [hep-ph/990649^ . 

[6] C. Bernard et al., Nucl. Phys. B (Proc. Suppl.) 4 (1988) 483. 

[7] E. Franco et al., Nucl. Phys. B317 (1989) 63; 
M. B. Gavela et al., Nucl. Phys. B306 (1988) 677. 

[8] S. R. Sharpe, A. Patel, Nucl. Phys. B417 (1994) 307; 

T. Bhattacharya, R. Gupta, S. R. Sharpe, Phys. Rev. D55 (1997) 4036. 

[9] M. Crisafulh et al., Phys. Lett. B369 (1996) 325; 
L. Conti et al., Phys. Lett. B421 (1998) 273. 

[10] JLQCD Collaboration (S. Aoki et al.), Phys. Rev. D60 (1999) 034511. 

[11] C. R. AUton et al., Phys. Lett. B453 (1999) 30. 

[12] A. Domni, V. Gimenez, L. Giusti, G. Martinelh, Phys. Lett. B470 (1999) 233. 

[13] M. Ciuchini et al., JHEP 9810 (1998) 008. 

[14] F. Parodi, P. Roudeau, A. Stocchi, Nuovo Cim. A112 (1999) 833; 



F. Caravaglios, F. Parodi, P. Roudeau, A. Stocchi, |hep-ph/0002171 



[15] M. Ciuchini, E. Franco, L. Giusti, V. Lubicz, G. Martinelh, Nucl. Phys. B573 (2000) 201. 

[16] L. Giusti, A. Romanino, A. Strumia, Nucl. Phys. B550 (1999) 3; 
R. Barbieri, L. J. Hah, A. Romanino, Nucl. Phys. B551 (1999) 93. 

[17] L. Maiani, G. Martinelh, G. C. Rossi, M. Testa, Nucl. Phys. B289 (1987) 505. 

[18] A. Donini, V. Gimenez, G. Martinelh, M. Talevi, A. Vladikas, Eur. Phys. J. CIO (1999) 121. 

[19] P. Hasenfratz, Nucl. Phys. B (Proc. Suppl.) 63A-C (1998) 53. 

[20] P. Hasenfratz, V. Laliena, F. Niedermayer, Phys. Lett. B427 (1998) 125. 

[21] H. Neuberger, Phys. Lett. B417 (1998) 141; 
H. Neuberger, Phys. Lett. B427 (1998) 353. 

[22] P. Hasenfratz, Nucl. Phys. B525 (1998) 401. 

[23] M. Liischer, Phys. Lett. B428 (1998) 342. 

[24] P. H. Ginsparg, K. G. Wilson, Phys. Rev. D25 (1982) 2649. 

22 



[25] T. Chiu, C. Wang, S. V. Zenkin, Phys. Lett. B438 (1998) 321. 

[26] D. Kaplan, Phys. Lett. B288 (1992) 342; 
Y. Shamir, Nucl. Phys. B406 (1993) 90; 
Y. Shamh^ |hep-lat/ 000302^ . 



[27] P. Hernandez, K. Jansen, M. Liischer, Nucl. Phys. B552 (1999) 363. 

[28] S. Capitani, M. Gockeler, R. Horsley, P. E. L. Rakow, G. Schierholz, 
Phys. Lett. B468 (1999) 150; 

S. Capitani, M. Gockeler, R. Horsley, P. E. L. Rakow, G. Schierholz, 
Nucl. Phys. B (Pr oc. Suppl.) 83-84 (2000) 893; 
S. Capitani et al., |hep-lat/000700^ . 

[29] H. Neuberger, Phys. Rev. Lett. 81 (1998) 4060; 

R. G. Edwards, U. M. Heller, R. Narayanan, Nucl. Ph ys. B540 (1999) 4 57; 
A. Bode, U. M. Heller, R. G. Edwards, R. Narayanan, |hep-lat/9912043| ; 



P. Hernandez, K. Jansen, L. Lellouch, hep-lat/OOOlOOS ; 



R. Narayanan, H. Neuberger, |hep-lat/0005004 



S. J. Dong, F. X. Lee, K. F. Liu, J. B. Zhang, |hep-lat/000600^ ; 



L. Giusti, C. Hoelbhng, C. Rebbi, Nucl. Phys. (Proc. Suppl.) 83-84 (2000) 896 and in prepa- 
ration. 

[30] M. Beneke, G. Buchalla, I. Dunietz, Phys. Rev. D54 (1996) 4419 and references therein. 

[31] C. Dawson et al., Nucl. Phys. B514 (1998) 313; 
M. Testa, JHEP 9804 (1998) 002. 

[32] M. Ciuchini, E. Franco, G. Martinelh, L. Reina, L. Silvestrini, Z. Phys. C68 (1995) 239. 

[33] G. Martinelh, C. Pittori, C.T. Sachrajda, M. Testa, A. Vladikas, Nucl. Phys. B445 (1995) 81. 

[34] C. Alexandrou, E. Follana, H. Panagopoulos, E. Vicari, bep-lat/0002010| . 



[35] S. Capitani et al., Nucl. Phys. B570 (2000) 393; 

S. Capitani et al., Nucl. Phys. B (Proc. Suppl.) 83 (2000) 232. 

[36] G. Martinelh, Phys. Lett. B141 (1984) 395. 

[37] M. Ciuchini et al., Nucl. Phys. B523 (1998) 501. 

A. J. Buras, M. Misiak, J. Urban, ^lep-ph/OOOSTSS . 



[39] M. Liischer, S. Sint, R. Sommer, H. Wittig, Nucl. Phys. B491 (1997) 344. 

[40] F. Niedermayer, Nucl. Phys. (Proc.Suppl.) 73 (1999) 105. 

[41] M. Ishibashi, Y. Kikukawa, T. Noguchi, A. Yamada, Nucl. Phys. B576 (2000) 501. 

23 



[42] S. Capitani, |hep-lat/0005008 . 

[43] V. Gimenez, L. Giusti, F. Rapuano, M. Talevi, Nucl. Phys. B540 (1999) 472. 

[44] A. Gonzalez- Arroyo, F. J. Yndurain, G. Martinelli, Phys. Lett. B117 (1982) 437; Erratum-ibid. 
B122 (1983) 486. 

[45] S. Capitani, M. Liischer, R. Sommer, H. Wittig, Nucl. Phys. B544 (1999) 669. 



[46] D. Becirevic at al., |hep-lat/00050T3 



24 



